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Abstract. The 71(1520) resonance has been recently studied in a unitarized coupled channel formalism 
with ttS(1385), KE(1530), KN and nS as constituents blocks. We provide a theoretical study of the 
predictions of this model in physical observables of the pp — » pK + K~p and pp — » pK + ir°n° A reactions. In 
particular, we show that the ratio between the tt°tv°A and K~p mass distributions can provide valuable 
information on the ratio of the couplings of the /1(1520) resonance to 7r 17(1385) and KN that the theory 
predicts. Calculations are done for energies which are accessible in an experimental facility like COSY at 
Jiilich or the developing CSR facility at Lanzhou. 
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1 Introduction 

, The low-lying negative parity resonances J P = 1/2- , 
3/2~, have recently attracted much attention as many 
of them can qualify as dynamically generated resonances 
from the interaction of mesons and baryons. In particular, 
much progress has been done interpreting the low lying 
1/2 - resonances as dynamically generated from the inter- 
action of the octet of pseudoscalar mesons with the octet 
of stable baryons |TT|2j|^]P||$||^j which allows one to make 1 
predictions for resonance formation in different reactions 
. One of the surprises on this issue was the realization 
that the .4.(1405) resonance is actually a superposition 
of two states, a wide one coupling mostly to irZJ and a 
narrow one coupling mostly to KN 0,8 9 . The perfor- 
mance of a recent experiment on the K~p — > ir°it°E 
reaction |10j and comparison with older ones, particu- 
larly the ir~p — > K°tt£ reaction has brought evi- 
dence on the existence of these two .4(1405) states |12| . 
Extension of these works to the interaction of the octet 
of pseudoscalar mesons with the decuplet of baryons has 
led to the conclusion that the low lying 3/2~ baryons 
are mostly dynamically generated objects [T3HT1| . One of 
these states is the yl(1520) resonance which was generated 
from the interaction of the coupled channels 7rZ'(1385) and 
KS(1530). From the experimental point of view this res- 
onance is of particular interest in searches of pentaquarks 
in photononuclear reactions |15U16| in jp — > K + K~p and 
7G? — > K + K~np. Since getting signals for pentaquarks 



involves cuts in the spectrum and subtraction of back- 
grounds, the understanding of the resonance properties 
and the strength of different reactions in the neighbor- 
hood of the peak becomes important in view of a correct 
interpretation of invariant mass spectra when cuts and 
background subtractions are made. From the theoretical 
point of view, in the studies of Refs. |13U14| . the .4(1520) 
is build up from the ^17(1385) and KS(153Q) and couples 
mostly to the first channel, to the point that, in this pic- 
ture, the state would qualify as a quasibound ir£* state. 
Indeed, the nominal mass of the vl(1520) is a few MeV be- 
low the average of the ir£* mass. However, the PDG [T7] 
gives a width of 15 MeV for the .4(1520), with branching 
ratios of 45% into KN and 43% into ttS, and only a small 
branching ratio of the order of 4% for ir£* which could be 
of the order of 10% according to some analysis JHj which 
claims that about 85% of the decay into 7T7ryl is actually 
ttS*. The association of irirA to ttS* in the peak of the 
.4(1520) is a non trivial test since one has no phase space 
for ttS* excitation and only the width of the S* allows 
for this decay, hence precluding the reconstruction of the 
S* resonant shape from the 7ryl decay product. 

One step forward in the understanding of the .4(1520) 
resonance has been possible thanks to a recent reaction 
K~p — > tt°tt°A, experimentally performed in Ref. fT!|| and 
theoretically studied in Refs. |20U21| . The reaction pro- 
ceeds mostly via K~p — > 7r°(i;* ) — ► 7r°(7r°/l). This is seen 
at energies of the K~ such that y/s > M 7r ox'.o, where the 
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reconstruction of the ir°A invariant mass produces the E*° 
resonance shape [THll^ . However, at energies of the K~ 
where y/~s ~ M_4( 152 o), the U*° is only produced through 
the tail of the resonance. Yet, a formalism using explic- 
itly the S* a propagator allows us to study properly the 
K~p — > tt°tt°A reaction assuming the same tt S*° as the 
dominant mechanism. 

On the other hand, the large branching ratios to KN 
and xS, of the order of 90% together, indicate that the 
KN and nE channels must play a relevant role in building 
up the 71(1520) resonance. In the work of Refs. |2()U21j 
this problem was tackled by performing a coupled channel 
analysis of the 71(1520) data with nE*, KS*, KN and wS, 
the first two channels interacting in s-wave and the last 
two channels in d-wave to match the 3/2 _ spin and parity 
of the 71(1520) resonance. In Ref. [23 it was also shown 
that although the ttU* remains with the largest coupling 
to 71(1520), its strength is reduced with respect to the 
simpler picture of only irZJ* building up the resonance, 
and at the same time there is a substantial coupling to 
KN and ttS which distorts the original irZJ* quasibound 
picture and makes the KN and ttS channels relevant in 
the interpretation of different reactions. 

In the present work we propose two reactions to test 
the non-trivial predictions of the unitarized coupled chan- 
nel model of Ref. [23 regarding the couplings of the 71(1520) 
resonance to the different channels. These reactions are 
PP ~ * pK + K~p and pp — > pK + Tr°ir°A close to the 71(1520) 
threshold. Particularly, we show that a measurement of 
the ratio between the 7r°7r°7l and K~p mass distributions 
for pp — > pK + 7r 7r°7l and pp — * pK + K~p reactions re- 
spectively is an excellent model independent test of the 
ratio between the 71(1520) coupling to 7r 17(1385) and to 
KN. 

The structure of the paper is as follows. In section [21 
the unitarized coupled channel model of Ref. [21] is sum- 
marized. In section[3thc dominant mechanism to produce 
the 71(1520) in p-induced reactions and the implementa- 
tion of the unitarized coupled channel model into it is 
described. Finally, in section0]we show the results for the 
mass distributions in the two reactions studied. 



2 Summary of the unitarized coupled channel 
formalism 



In Ref. |23 the 71(1520) resonance was studied within a 
coupled channel formalism including the ttS*, KS* in s- 
wave and the KN and ttS in d- waves. Unitarity was im- 
plemented by means of the Bethe-Salpeter (BS) equation 
in the evaluation of the different scattering amplitudes, 
which reads 



T = V + VGT => T = [1 - VG^V. 



Table 1. Loop subtraction constants and parameters of the 
potentials. The 713 and 714 are given in units of GeV -3 ; and 
733, 744 and 734 in units of GeV -5 . 



do 


«2 713 


714 


733 


744 


734 


-1.8 


-8.1 98 


110 


-1730 


-730 


-1108 



The kernel of the BS equation is given by the following 
transition potentials 21 



V = 



Cxi{k° x 
C 2 i(fc 2 ° 



k?) Cu(fc? 
*?) C 22 (k° 



fc 2 °) 
fc 2 °) 



713 93 



713 93 

714 94 



7339! 

734 93 94 



714 94 



734 93 94 

74494 



where the elements 1, 2, 3 and 4 denote irZJ* 
and ttS channels respectively. In Eq. J2J), 



KS* 



2Vi 



y/[s - (Mi + m,) 2 ][s - (Mi - to,) 2 



(2) 
KN 



(3) 



k? 



s-M?+m? 



and Mi(nii 



The coefficients G\j are Cn 



is the baryon(meson) mass. 

_ Ve 



J- 



C: 



C21 =C 12 = ^ and 
where / is 1.15/ ff , with f n (= 93 MeV) the 
pion decay constant, which is an average between f„ and 
fx as was used in Ref. [2] in the related problem of the 
dynamical generation of the 71(1405). In Eq. QJ G stands 
for a diagonal matrix containing the loop functions in- 
volving a baryon and a meson which are regularized by 
means of two subtraction constants [2J|: one for the s- 
wave channels (do) and another one for the d-wave chan- 
nels (a 2 ). The treatment of Eqs. (1) in Refs. PUII^T] relies 
upon a dispersion relation on T^ 1 which allows the on- 
shell factorization out of the loops of the kernel, V, of the 
Bethe-Salpeter equation [21122] • The matrix elements Vi%, 
V12, V21 and V22 come from the lowest order chiral La- 
grangian involving the decuplet of baryons and the octet 
of pseudoscalar mesons, as discussed in Ref. |23II13II14| . 
The unknown parameters in the V— matrix, as well as the 
subtraction constants of the loop functions, were obtained 
by a fit to KN — > KN and KN — > irU partial wave am- 
plitudes. The values obtained are shown in Tabled 21 . 
Despite the apparent large number of parameters, it is 
worth mentioning that the dominant potentials are Vn, 
V12, V21 and I7 22 which have no freedom. Let us recall that 
we do not include ttS* amplitudes in our fit. Hence, the 
amplitudes involving this channel are genuine predictions 
of the theory and are shown in Fig. ^ 

From the scattering amplitudes we can also extract the 
effective couplings of the 71(1520) to all the different chan- 
nels, which can be obtained via an analytic continuation 
of the amplitude in the complex plane. Near the pole we 
may write up to regular terms 



- M A{1520) + 2-T4(1520)/2 



(4) 
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nL*~7lE* JlE^KN 7tE*-*Jl£ 
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Fig. 1. Unitary amplitudes involving the nE* channel. From 
left to right: irE* — » -nE* , irE* — > K N and ttE* — > nE. 

Table 2. Couplings of the yl(1520) resonance to the different 
channels 



9i 


92 


93 


54 


0.91 


-0.29 


-0.54 


-0.45 



from where we have 

-^1(1520) \T % j{M A ( lb2 Q))\ 2 . . 

where M^ 152 o) is the position of the peak in \Tij\ 2 and 
^1(1520) = 15.6 MeV. The couplings obtained are shown 
in Tabled 



3 Proton-induced /1(1520) production 

The amplitudes involving TtE* channels can be investi- 
gated in particular in those reactions where this chan- 
nel plays a significant role. In Ref. [21], these amplitudes 
were checked in the K~p — > Ann, jp — > K + K~p, jp — > 
K + ir°ir°A and ir~p — > K°K~p reactions, leading to a 
good reproduction to the experimental results. 

In the present work we apply the model to the pp — > 
pA'+ (4(1520)) -> pK+(tt tt A) andpp -► pjf+ (4(1520)) - 
pK + (K~p) reactions at energies slightly above the 4(1520) 
production threshold, which are attainable at a facility like 

cosy [21ES1. 

As explained in Appendix^] only a single partial wave, 
with the initial pp pair in the 1 D 2 1 , contributes to the pro- 
cess pp — ► K + pA(1520) near threshold under the assump- 
tion that the 4(1520)KN system in the intermediate state 

1 Here we use the standard notation for the NN partial 
waves: 2S+1 Lj, with S, L, J for the total spin, the angular 
momentum and the total angular momentum. 



is in an s-wave in all subsystems. In addition, since the 
production is characterized by a large momentum trans- 
fer, we can assume the production operator to be (largely) 
independent of the final momenta, that are constrained to 
small values because of the chosen kinematics. Thus the 
only significant source of an energy or momentum depen- 
dence with respect to the final particles should be their 
various final state interactions. Therefore, to deduce infor- 
mation on the final state interactions in large momentum 
transfer reactions, no detailed knowledge of the produc- 
tion operator is necessary (for a recent review on these 
issues we refer to Ref. and the references therein). 

There are various ways to derive the relevant matrix 
elements. The one most commonly used is based on the 
irreducible tensor techniques and the corresponding for- 
mulas are given in Appendix B. Here on the main text, 
on the other hand, we will describe a method for the con- 
struction of matrix elements that is more transparent and 
allows us to easier explain the relevant physics. Needless 
to say that the final results in the two approaches are 
identical. 

According to the selection rules given above the tran- 
sition pp — > 4(1520)Ap can be parameterized by a single 
constant C and a fixed operator structure which reads for 
non-relativistic final states 

W = C ((m^ (1520) • p)(a ■ p)(T 2 u*^j 4> K (uJct 2 u p ) , (6) 

where it a and u p denote the spinors for the .4(1520) as well 
as the nucleons. Note, by construction er-it 71(1520) vanishes. 
It is this identity that ensures that it is only the pp D wave 
that contributes to the above transition operator. For the 
initial momentum we use p, and er denotes the standard 
three vector of Pauli matrices. In the expressions below 
we will omit the spinors from the initial state to simplify 
notations. 

To come to the corresponding expressions for the full 
transitions — including the decay of the 4(1520) — we need 
to contract the transition operator W with further oper- 
ators, namely the spin transition operator S normalized 

as 

SiSj = - (2S lj - ie ljk o k ) , 

as it occurs in the non-relativistic version of the Rarita- 
Schwinger propagator, as well as the relevant vertices and 
propagators. As a result for the decay matrix element for 
the chain 

pp -► K + pA{imQ) -> K + p(ir°£*°) -> K + p(Tr°[ATr }) 
we need to replace UM 152 o)i m Eq. © by 

= -± gi l^-u A (p 3 )(S-p 1 )SlG A ,(s A *)G s ,(s 1 z*) , 

V O m ?T 

where we used for the transition 4(1520) — * £*ir a con- 
stant matrix element and for S* — > Att the standard 
vertex S ■ p\. In the expression — gi/vo is the effective 
coupling for 4* — > ir°Z!* a including the isospin factor. 
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Here Ga* denotes the propagator of the ^4(1520), im- 
plicit in the dashed circle in Fig. |3 which is a function 
of sa- = (pi +P2 +Pz) 2 - Gz* is the propagator of the S* 
being a function of s^.* = {p 2 +P3) 2 , where we use the 
standard Flatte-parametrization 

G E .(s E .) = — 



My 



Note, to make contact with the T-matrices derived in the 
previous section we need to replace giGyi* by the corre- 
sponding channel T matrix. We come back to this below. 

After some algebra we get for the complete matrix el- 
ement in the irirA channel 

V 3 m 7T 

x {V(pi) i G E .(8 1 E .)+Vfa) i G E .(sl.)} ,(7) 

where V{k)i — (k ■ p)pi — \p 2 ki projects on the initial 
D-wave, as already explained above. One observes that 
the Ap system in the final state occurs solely in the spin 
triplet in line with the findings of Appendix B. 
On the other hand, for the reaction chain 



PP 



K+pA(l520) -> K+p(K-p) 



we need to replace tt^n 52 o) j m Eq. © by 

r K = ■ p 3 )(S ■ p 3 )St . (8) 

We get after some standard manipulations 

A Kp (Pi,P2) = -^ C 93 {u(p 1 )a 2 u(p 2 )(V( P3 ) - P3 ) 

x (Ga*(si) + G A '(s2)) 

+iu(p 1 )<T k <T2u{p 2 )(P3 X p)fe(f> 3 • p) 

x {G a »{sx)-Ga»{s 2 ))}{$) 

Here the first term contains the two proton system in a 
spin singlet state (e.g. 1 5'o), whereas the second term con- 
tains the pp system in a spin triplet state. The second 
term implies L = odd for the two protons. However, as ex- 
plained in Appendix IbI the energies of the proton are such 
that only L — is relevant. In addition, Ga* (si)-Ga* (s 2 ) 
is much smaller than the sum. Therefore we will only in- 
clude the first term in what follows. 

As long as we study ratios only, the only unknown in 
the above expressions, C, drops out and all distributions 
turn out to be predictions of the model. However, to make 
contact between the coupled channel T-matrices of the 
previous section and the expressions just derived we use 
a particular production mechanism as shown in Figs. 
and |3| Thus the transition operator parameterized as C 
contains — in this model — the KKNN vertex, where the 
outgoing K + is produced, as well as a Kp vertex on the 
other nucleon. We may thus pull the latter out of the 
definition of C and write 

-±fm = c. 




Fig. 2. Mechanism considered for pp — > pK + (yl(1520)) 
P K + (tt°ty°A) 

P P 

/\P 2 APl 




KT 




Aq 2 
(-m) 



Fig. 3. Mechanism considered for pp — ► pK + (A(1520)) 
pK+{K-p) 



With this definition it is straight forward to show that 



C 



1 

W 



(p2-*°) 



1 



k 2 - 



(10) 



for the transition operator and 



(5l) G A ( 15 2 ) 



( 9-ap" 



\V3 
( 93P 2 



( 93Pl \ r 

W§7 Ga(1520) V Vs 



— N — >tt£* 7 (11) 

— ^kn^kn (12) 



for the meson-baryon T-matrices. This completes the eval- 
uation of the production matrix elements. 

By the described method the antisymmetrization of 
the two proton states has already been taken into account 
since only allowed partial wave were considered. The sym- 
metrization of the two identical ir° has to be considered 
by adding the contribution obtained by changing p\ «-> p 2 
(see details in Appendix [HJl . Hence, |A™p| 2 to be used in 
the cross section (c.f. Eq. Ijl9|l in the appendix) is given 

by 



\Amp\ 2 



\A™ A ( Pl ,p 2 )+A™ A (p 2 , Pl )\' . (13) 



For the calculation of the cross sections one has to evaluate 
a four- and five-body phase space. The method we used 
is given in detail in Appendix lAl 
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Since we have two baryons in the final state which are 
different for each reaction, the different modification of the 
cross section due to the final state interaction (FSI) of the 
AN or pp subsystem can be relevant. We have taken this 
into account by means of the ordinary factor, equivalent 
to the inverse of the Jost function 2 . For the Ap interaction 
we use [2HJ for the factor that multiplies |T| 2 in the phase 
space 



\C. 



FSI\ 



q 2 +l3 2 



(14) 



with 



a = (1 - y/1 - 2r/a)/r, = (1 + y/l - 2r/a)/r (15) 

and q the A momentum in the pA rest frame. In Eq. (|15|) , 
a and r are the scattering length and effective range of 
the s— wave Ap amplitude in S = 1 which is the state we 
have. In the calculations we use a — (— 1.4 ± 0.5) fm and 
r = (4 ± 1) fm and we shall estimate the uncertainties 
induced from these errors. 

For the pp final state interaction in the pp — » pK + K~p 
reaction we use the same expression but with a = —7.8 fm 
(which already accounts for interference with the Coulomb 
force [HO]) and r = 2.79 fm with significant smaller errors 
than in the Ap case. This simple prescription is accurate 
to better than 10% in the region of energies relevant here 
|31| compared to more elaborate formulas as, e.g., the one 
given in Ref. t 32j and is therefore sufficiently accurate for 
our purposes. 



4 Results 

First we show in Fig. [2] the results for the ir°ir°A invari- 
ant mass distribution of the pp — > pK + ir°ir°A reaction 
and the K~p mass distribution for the pp — > pK + K~p 
(dashed and solid lines respectively) for an incident pro- 
ton momentum of 3.7 GeV which is the highest available 
at COSY. This is also the nominal energy for protons of 
the developing CSR facility at Lanzhou (China). The dif- 
ferent shapes below the peak for the two reactions can be 
understood from the presence of the S* propagator in the 
pp — ► pK + tt q tt A reaction, since only its quickly decreas- 
ing tail enters in the evaluation of the matrix element. 

Although the production rates with this model are of 
the order of experimental cross sections measured in the 
pp — > ppK + K~~ reaction |25|. other mechanisms can be at 
play. In addition, we did not include the initial state inter- 
action that is expected to lead to a sizable reduction of the 
cross section [33], but would not modify the ratio of cross 
sections. Because of this we focus on the shape and rela- 
tive magnitude of our results that are model independent 
as explained above. In an analogous way to what was stud- 
ied in Ref. |2J for the photoproduction case, we propose 
to measure the ratio between the 7r 7r°yl and K~ p mass 



2 Although the Jost function should not be used for the 
extraction of scattering parameters from production reac- 
tions 27 , it still gives a for our case sufficiently accurate 
parametrization for the effect of the final state interaction |28) . 




1480 1500 
M (MeV) 



4. Solid line: K p invariant mass distribution for the 



Fig. 

pp — > pK + K~p reaction. Dashed line: ty v ty v A invariant mass 
distribution of the pp — > pK + ir°-K A reaction. Incident proton 
momentum: 3.7 GeV. 



distributions of the pp 
reactions respectively 



• pK + n a Tr A and pp — > pK + K~ 



R = 



da 



UU pp—>pK+ 7T 1 



da 



pp- 



(16) 



->pK+K-p/ dM K - p 
Up to phase space, FSI and known numerical factors, 
the ratio R is proportional to (\T- K £*->- k E'\/\Tk£*^kn\) 2 
which, at the yl(1520) peak position, is (31/53) 2 = 2.8 (see 
Eq. (@J). With the full model, the value obtained for R at 
the peak position is 



R~ 0.5 ±0.2. 



(17) 



The uncertainty given comes from the errors of a and r 
in Eq. I|15|) when considering the FSI of the final baryons, 
which is the largest source of uncertainty in our model. 
Had we not considered the FSI we would have obtained 
R ~ 0.9. Actually the effect of the FSI is to increase by 
a factor of ~ 3.6 the K~p mass distribution and a fac- 
tor ~ 2 the tt°w°A mass distribution at the peak in the 
pp — > pK + K~p and pp —* pK + ir a ir°A reactions respec- 
tively. These factors emerging from the final state inter- 
actions are strongly dependent on the beam energy. 

An experimental measure of this ratio R would provide 
a good test of the unitarized coupled channel model since 
the value of g\ and g$ are non-trivial genuine predictions 
of this theory. Note that the close connection between the 
nature of a resonance and its effective couplings to the 
decay channels was derived already in Ref. 34 for bound 
states and extended to inelastic resonances in Ref. 35 . 



5 Summary 

We have studied some of the consequences and predic- 
tions of a unitarized coupled channel approach to the 
71(1520) resonance through the pp — ► pK + K~p and pp — > 
pK + ir°ir°A reactions. The model relies upon a coupled 
channel formalism implementing unitarity through the Bcthc- 
Salpeter equation by means of which the resonance struc- 
ture appears naturally. One of the genuine predictions of 
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this theory are the values of the effective couplings of the 
.4(1520) resonance to the relevant channels. We propose 
that the ratio between the tt tt A and K~p mass distribu- 
tions for the pp — ► pK + ir°ir°A and pp — > pK + K~p reac- 
tions respectively, can provide a test of the ratio between 
the coupling of the 4(1520) resonance to the 7ri7(1385) 
and to KN channels. Such a ratio would help to un- 
ravel the nature of the 71(1520) and its coupling to meson- 
baryon components. 



is Lorentz invariant, hence we can evaluate it in any de- 
sired frame. By convenience we choose a frame where 
P3 +P4 = 0. In this frame the 5 function of energy conser- 
vation is simple since M| 4 = (p3 + pa) 2 = (qi + q 2 — pi — 
P2~Ph) 2 = {q° + q% - Pi - P2 - P%) 2 and hence + 



P1-P2-W3-U4-P5) = S(M 3i ~ y/pl + m§ - y/pl + ml) 
is only a function of |p 4 | and hence can be use to evaluate 
the \p±\ integral, which gives 
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1 \p'A 



M: 



34 



A Five- and four-body phase space 

Here we explain in detail the steps followed to simplify the 
integrals of the phase space in the evaluation of the cross 
section for the 

P(gi)p(ffi0 -> p(p 5 )K+(p 4 )ir ( Pl )ir (p 2 )A(p 3 ) (18) 

reaction (see Fig. |21 for the detailed definition of the mo- 
menta) with masses Mi, M%, 7715, 7714, mi, mi and m 3 
respectively. 

The cross section for the process is given by: 



2SM 1 M 2 



y/s 2 - 2s{M 2 + M 2 ) + (M 2 - M 2 ) 2 

d 3 pi f d 3 p 2 f d 3 p 3 f d 3 pi f d 3 p 5 



(2tt) 3 J (2tt) 3 J (2tt) 3 
1 1 m 3 1 TO5 
2u>i 2u>2 0J3 2oj 4 UJ5 

<!:l ' J^Pi) A "' P ' 



(2tt) 3 J (2tt) 



x (27r)V( gi 1 X 



(19) 



where 5 = 1/2 for the symmetry of 7r°7r° andpf = (wi,pi). 

The evaluation of the S of energy conservation for so 
many particles can be a difficult task. However it can be 
extremely simplified by using the following procedure. 

The quantity 



A 



d 3 pi f d 3 p 



1 



(2tt) 3 J (2tt) 3 lu 3 lo 4 

X {2*f5\q 1 +q 2 -Y J Pi)\ Am Pf 
d 3 p 4 1 



(27r) 3 

x (2n)6(q° 1 +q2-J2o Jl )\Amp\ 2 



(20) 



I Amp 1 2 (M34 — m 3 — 777,4 ) 
(21) 

where \p' 4 \ = \ 1 / 2 (M 34 ,m 3 ,m 2 : )/2M 3 4 is the momentum 
in the P3 + J>4 = frame. 

The final expression for the cross section is 

S f Af 1 M 2 TO3m5 
4(2 7 r) li y^ 2s{M 2 + M 2 ) + (M 2 - M^f 

x J d 3 p x J d 3 p 2 J d 3 p 5 J dcosd' 4 J d<j)' 4 



W1W2W5M34 



6*(M34 — 777,3 — mn)\AmpY 



(22) 



When evaluating the amplitude all the momenta have to 
be evaluated in the overall CM frame, hence one has to 
boost back the p' 4 momentum to this system with the 
following expression: 



P4 =P 4 



/ p0 



- 1 



Pi 



' P 



Jo 
Pi 



VM34 7 \p\ 2 ' M34 

with P = — pi P2 — Ps and P° = yfs — p\ — p\ — p% 
The four-body phase space for the 



(23) 



p(qi)p(q2) -> p{p2)K + {p 4 )K (p 3 )p(pi 



(24) 



reaction can be evaluated in an analogous way to the pre- 
vious one and gives 

SM\M2m\m2 
2(2i:) s y^ 2 ^2s(M 2 + M 2 ) + (M'( - M 2 f 



d 3 p! J d 3 p 2 J dcos 6» 4 J d<j)' 4 
x ©(M34 — m 3 — mi)\Amp\ 2 , 



wiw 2 M 34 



(25) 



where 5=1/2 since the final protons are identical parti- 
cles. The boost of p' 4 to the overall CM frame is done with 
Eq. J23 but now P = p ± - p 2 and P° = ^fs - p\ - p\. 



B Alternative derivation of the amplitudes 

Here we explain in detail the analytic evaluation of the 
amplitudes for the pp — > pK + 7r°7r°yi and pp — > pK + K~p 
reactions within the method of irreducible tensors, which 
is more familiar to many readers than the technique de- 
scribed in the main text. 
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Let us start by determining the allowed quantum num- 
bers of the initial and final states, which are independent 
of the internal mechanisms, which will be of help in simpli- 
fying the evaluation of the amplitudes later on. Let us con- 
sider the general process pp — > K + pA(152Q). Since the re- 
action is calculated at energies close to threshold, the only 
final total angular momentum allowed is L — 0. Hence the 
final J p possibilities are 1 + or 2 + , since J p (K + ) = CP, 
J p {p) = 1/2+ and J p (/l(1520)) = 3/2-. The spin of the 
initial pp pair can be or 1 while L has to be even for 
parity reasons. But, since the initial protons are identical 
fermions, L + S + 1 has to be odd. And, since, 1=1, then 
L + S has to be even. Hence the only possibility to match 
final 1 + or 2 + is that the initial protons are in L = 2 and 
S = 0. In summary, the initial pp state has the following 
quantum numbers: L = 2, S = 0, P = +, independent of 
the internal dynamics and the 71(1520) decay products. 

Let us consider now the quantum numbers of the pA 
system in the pp -> pK+(A(1520)) -> pK + (ir Q £*) -> 
pK + (ir°ir°A) reaction. Considering the parity of the final 
particles and the fact that the 71(1520) decay into irU* is 
in s-wave and the decay of the £* into irA is in p-wave, 
the angular momentum of the pA system has to be even to 
match the global parity +. Given the typical momenta of 
the final A at the energies of concern in the present work, 
the only possibility is L(pA) = 0. On the other hand the 
spin of the pA system has to be S(pA) = 1 in order to give 
total J — 2 when combined to the p-wave of the S* — » ttA 
decay in line with the formalism as described in the main 
text. 

Regarding the pp — > pK + (71(1520)) — » pK + (K~p) re- 
action, the kinematics of the pp final state is that the pro- 
ton coming from the 71(1520) decay has about 240 MeV/c 
momenta since the 71(1520) is produced essentially at rest, 
and the other proton is also basically at rest in the to- 
tal CM frame. This gives the protons a kinetic energy of 
about 8 MeV in their CM frame where only relative S- 
waves are relevant. Hence, the final pp system only can 
be in L{pp) = at the energies considered. On the other 
hand the spin is S(pp) — since L + S + I has to be odd. 

Let us evaluate the amplitude for the pp — > pK + ir°ir°A 
reaction of Fig. [5] The amplitude for the vertex pK + K~p 
is given by 



it 



1 



(26) 



The unitarized KN — > ttS* transition has the following 
form 20^ 



itK-p-+TT Q £'° — —i 



WW C (2 2 2>- m > M 



x y 2 ,- m -M(fc)(-ir + "V47r, 

where C are Clebsch-Gordan coefficients. For the decay of 
the U* into tt°A, the vertex is 36 

-it = - f -^{\m"\S-p*\~M). (28) 



where fs'TrA/m^ — 9.61 x 10 3 MeV 1 and S is the total 
spin 3/2 to 1/2 transition operator given by: 



Sp = 







fp, 



■Pz 







Px-Wy 



(29) 



V2 



t = i 



Hence, the full amplitude is given by 

1 , n .n, 1 1 



2f 



k 2 - 



/s^ - M s , + ir E . ( v / ^ 7 )/2 



^ fs*nA 1 rp 
X —1 ; 



1 



a v^ttS* C{- 2 -;-m,M + m) 



x Y 2 ^ M (k)(^l) M+m V^(\m''\S-p 2 \^M) 

X S m ^ m SM,-m, (30) 

where s s * = (jp2 + Pa) 2 ■ 

Since the momentum of the exchanged kaon is approx- 
imately equal to the incoming proton momentum, we can 
do the following approximation: 



(31) 



On the other hand, the spin of the initial pp system 
is S(pp) — (antisymmetric) and the spin of the final 
pA system is S(pA) = 1 (symmetric), as discussed above. 
Therefore, the initial spin wave function has to be anti- 
symmetrized and the final pA spin wave function has to 
be symmetrized. Let us consider only the spin dependent 
part of Eq. I|30|l . The antisymmetrization of the initial spin 



gives 



1 

V2 



C(^2^;-io)<im"|5.p 2 |-i)«5 m ,, +1/2 
,1 _ 3 



^2 2' 2 "2 1 r '2 



J_ /I 
V2 \ 5 



l -Q){\ m "\S-p^- + \)5. 

)S m ',+l/2 



m',-1/2 



Lrn"\S-p 2 \- 



(I m "|S-p 2 || + i)5 ro , ! _ 1/2 



(32) 



Now we have to symmetrize the latter expression for 
the final pA spin. For S z (pA) — gives 

~ ,3 1. 



1 1 /2 



5-5|5-pa|r-r)+0 



0+(^ 



2> 



+ m) 
(27) 



2 12 
5 V 3 P22 



(33) 



for S z (pA) = -1: 







V2V5 

= ~7E7E iP2x + ip2y) ' 



2' 



(34) 
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and for S z (pA) - 

V2V5 

1 1 



[~ 2 + - 2 \s- P2 \- 2 -- 2 ) 

■(P2x - iP2y)- 







V5 V6 

Hence, the final expression for the amplitude is 



(35) 



t = 



1 



3 2/2 



fc 2 - 



is* * j 



{ V2p' 2z S z = 



\iP2x ~ W'ly) S z = +1 
, -Wlx +W2y) S z = -1 



(36) 



and the symmetrization of the identical pions has to be 
considering by adding, for each S z , the amplitude chang- 
ing pi <-► p 2 . 

For the evaluation of the pp — > pK + K~p amplitude 
we need the unitarized K~p — > K~p amplitude which is 
given by 



, 2 3 

itn-p^K-p = -i^T RN ^ RN ^C{- 2 -;-m,M + m) 

M 



x F 2 ,_ TO _ M (fc)C(i 2 |m",M-m") 



x ^>V-m(P3)(-1) 



m. +m 



47T. 



(37) 



The rest of the procedure is analogous to the previous 
case but now the final pp spin is S(pp) = 0, as discussed 
above, which is antisymmetric and hence we have now to 
antisymmetrize the final pp spin wave function. 

The final expression of the amplitude is given by 

t{pi,P2) = 2^(P4- fc °) fc2 ^ m 2 7^(3 COS 2 °i~ l ) T KN^KN 

(38) 

where the symmetry of the final pp space wave function 
has to be considered by adding to the amplitude the same 
expression changing pi <-> p 2 . 
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